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Abstract 

The mixed-isospin vector current correlator, (0|T(V^fV^)|0) is evaluated us- 
ing both QCD sum rules and Chiral Perturbation Theory (ChPT) to one-loop 
order. The sum rule treatment is a modification of previous analyses necessi- 
tated by the observation that those analyses produce forms of the correlator 
that fail to be dominated, near q 2 = 0, by the most nearby singularities. 
Inclusion of contributions associated with the 4> meson rectify this problem. 
The resulting sum rule fit provides evidence for a significant direct to — > tttt 
coupling contribution in e + e~ — » 7r + 7r~. It is also pointed out that results 
for the g 2 -dependence of the correlator cannot be used to provide information 
about the (off-shell) ^-dependence of the off-diagonal element of the vector 
meson propagator unless a very specific choice of interpolating fields for the 
vector mesons is made. The results for the value of the correlator near q 2 = 
in ChPT are shown to be more than an order of magnitude smaller than those 
extracted from the sum rule analysis and the reasons why this suggests slow 
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convergence of the chiral series for the correlator given. 
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I. INTRODUCTION 



Non-electromagnetic isosopin breaking is well-established in many strongly interacting 
systems (e.g., splittings in the hadron spectrum, binding energy differences in mirror nuclei, 
asymmetries in polarized np scattering, binding energies and level splittings of light A hy- 
pernuclei [TTT] ) . In few-body systems, an important source of this breaking has been thought 
to be the mixing of isoscalar and isovector mesons appearing in meson exchange diagrams. 
In particular, the bulk of the non-Coulombic contributions to the charge symmetry breaking 
nn-pp scattering length difference and to the A=3 binding energy difference, and of the np 
asymmetry at 183 MeV, can be explained ||,[3[ using the value of p — uj mixing extracted 
from an analysis of e + e~ — > 7r + 7r~ in the p — uj interference region f4j|| . The plausibility of 
this explanation (which employs the observed mixing, measured at q 2 = m 2 , unchanged in 
the spacelike region q 2 < 0) has, however, recently been called into question by Goldman, 
Henderson and Thomas || who pointed out that, in the context of a particular model, the 
relevant p — uj mixing matrix element has significant q 2 - dependence. Subsequently, various 
authors, employing various computational and/or model framewords, have showed that the 
presence of such g 2 -dependence appears to be a common feature of isospin-breaking in both 
meson-propagator- and current-correlator matrix elements PHT6 . 



In the present paper we will concentrate on the isospin-breaking vector current correlator 
IV(g 2 ) = i [d 4 x e«*(Q\T(V>{x) C(0))|0) , (1.1) 



where 



V£ = (ujpu - d lfl d)/2 , V£ = (ujpu + d^d)/6 . (1.2) 
This correlator was first analyzed using QCD sum rules in Ref. [0 , and the analysis 



updated by the authors of Ref. |L2] who, in particular, stressed the g 2 -dependence of the 



correlator implicit in the results of this analysis. As will be shown below, a worrisome feature 
of the resulting fit is that the phenomenological representation of the correlator near q 2 = 
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is not dominated by the most nearby singularities, suggesting that some ingredient may be 
missing from the form chosen for this representation. This missing ingredient is identified 
below and it is shown that a reanalysis of the correlator, which includes it, rectifies the 
problem. The resulting correlator still displays a very strong (^-dependence, and, in addition, 
provides evidence for the presence of significant direct u — > nn coupling in e + e~ — > 7r + 7r~ . 
The behavior of the resulting correlator near q 2 = is then compared with that obtained 
from ChPT to one-loop. The latter is found to be more than an order of magnitude smaller 
than the former, the reason why this suggests the likelihood of a slow convergence of the 
chiral series for the correlator explained. 

The paper is organized as follows. In Section II, those features of the behavior of quantum 
field theories under field redefinitions relevant to attempts in the literature to relate meson 
propagators and current correlators are discussed, and it is explained why the freedom of field 
redefinition implies that (1) one cannot obtain off-shell information about the off-diagonal 
element of the vector meson propagator from the off-diagonal element of the vector current 
correlator without making specific choices for the vector meson interpolating fields, and (2) 
if one writes the off-diagonal element of the vector meson propagator as 

6 puJ {q 2 ) cannot, in general, be g 2 -independent. In Section III we return to the QCD sum rule 
analysis of the vector current correlator, first explaining why certain features of the existing 
analyses suggest the need for a modified analysis, and then performing this analysis. The 
results both correct the apparently unphysical features of the previous analyses and provide 
evidence for non-negligible direct uj — > 7r + 7r~ contributions to e + e _ — > n + n~ in the p — uj 
interference region. In Section IV, the correlator is computed to one-loop in ChPT, and 
the results compared to those obtained from the sum rule analysis. Implications for the 
discrepancy between the results of the two approaches are discussed there. Finally, in 
Section V, a brief summary of the main results of the paper is given. 
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II. CONSEQUENCES OF THE FREEDOM OF FIELD REDEFINITION 



Let us begin by clarifying the relation (or lack thereof) between the vector-meson- 
propagator and vector-current-correlator matrices. The former is an, in general, off-shell 
Green function, which we may think of as being associated with some low-energy effective 
Lagrangian, £ e //, in which the vector meson degrees of freedom have been made explicit. 
As is well-known [T^-2D| , the form of such a Lagrangian is not unique: if <fi and x are t wo 
possible field choices describing a given particle, related by = X^(x)j with -F(O) = 1, then 
C e ff[4>] and £' e ff[x] = ^efflx-^ix)] produce exactly the same experimental observables |Tj| . 
However, while the S-matrix elements of the two theories are identical, this is not true of the 
general off-shell Green functions. One is free to make field redefinitions of the form above 
(as is done, e.g., in order to obtain the canonical form of the effective Lagrangian for ChPT 
IT9H2T1 ) without changing the physical consequences of the theory; the Green functions, 
however, are not in general invariant under such field redefinitions. Useful pedagogical illus- 
trations of this general principle, for pion Compton scattering and the linear cr-model, are 
given in Ref. [ 22[ and Chapter IV of Ref. [^3| , respectively. In the case of interest to us, what 
this means is that, when we make a redefinition of the p, uj fields in C e ff[p, u], we generate a 
new effective Lagrangian, C e ^[p' , uj'], the Green functions of which are, in general, different 
from those of £ e // (though when we piece such Green functions together to form S-matrix 
elements, these differences produce no net effect). The off-shell behavior of the vector meson 
propagator is thus dependent on the particular choice of fields used to represent the vector 
mesons (the choice of "interpolating field"). It is not a physical observable. In contrast, the 
vector current correlators nj(g 2 ) = i J d 4 x e iq - x (0\T(V°(x) V; 6 (0))|0) are, in fact, physical 
objects, independent of interpolating field choice. The spectral functions for LT^ and ITj^, 
are, for example, accessible from a combination of r~ — > u t tt^tt° and e + e _ — > tttt , tttttttt 
data, and that for IL^, could in principle be obtained from a careful analysis of the deviation 
of the ratio of the differential decay rates for r _ — > u T 7i^7T and e + e~ — > 7r + 7r~ from that 
predicted by isospin symmetry. As such there can be no general (i.e. valid for all choices of 
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interpolating field) relation between the correlator and propagator matrices. This point is 



the source of some confusion in Ref. |T2[ where an attempt is made to obtain the off-shell 
propagator based on an analysis of the correlator. 

Before proceeding to the reanalysis of the correlator, let us be more precise about the 
problems with the interpretation of the results of Ref. |T2]], in the light of the above com- 
ments. The authors begin by writing a general form for the spectral function of the corre- 
lator: 

ImIV(g 2 ) = A lmU^(q 2 ) + A, ImlLff + • • • (2.1) 

where the superscripts on the RHS should, for the moment, be taken only as labelling the 
region of spectral strength, and where + • ■ • refers to all other contributions (we return to 
this below). Eqn. (2.1) is, of course, completely general. The authors of Ref. [IT2J, however, 



then identify A Q with rr?m^l 9p9u) where g PtW are the vector meson decay constants, defined 
by 

9 

Tfl 

(0\V^\p,co) = -f^e,, (2.2) 
9p,u> 

and 11^ with the off-diagonal element of the vector meson propagator. This amounts 
to assuming that the isospin-unmixed I = 1 p state, p^°\ couples only to V? and the 
isospin-unmixed I = uj state, only to V p , the isospin-breaking contribution to IL^ of 
Eqn. (2.2) from the p,u region then resulting solely from the p^-u^ ' mixing in the meson 
propagator. In this interpretation, fixing the imaginary part of the correlator in the p, u 
region (via the sum rule analysis) allows one to obtain the isospin-breaking parameters of the 
imaginary part of the vector meson propagator, and, via a dispersion relation, the behavior 
of the off-diagonal element of the propagator off-shell. However, as explained above, such 
a possibility is excluded on general grounds. The problem with going from Aq and IT^ of 



Eqn. (2.2) to the interpretation of these quantities in Ref. [12] is that, not one, but three 
sources of isospin breaking exist in the contributions to Tl pu from the p, uo region: that due 
to p^-uj^ mixing (discussed above), that due to the direct coupling of V£ to cu^°\ and that 
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due to the direct coupling of Vjf to p^> . The same AI = 1 strong operator which gives rise 
to non-zero p^-tu^ mixing will also necessarily give rise to the latter two couplings. These 
couplings would be described by new isospin breaking parameters, 0^ and <f)^ p , 

(0|y;| p (o)) = 

Qui 

2 

(0|F» = -^>% (2.3) 
9p 

where <J)^ P , (/>^ u are, in general, g 2 -dependent, and also interpolating-field-dependent 
off-shell. Thus, off-shell, the p-u region contribution to IT^ depends not only on the 
(interpolating-field-choice-dependent) isospin-breaking parameters of the off-diagonal ele- 
ment of the vector meson propagator, but also on the (interpolating-field-choice-dependent) 
isospin-breaking parameters <p^ p , (j)^^ . The total contribution is independent of the in- 
terpolating field choice, but the individual contributions are not. One is, of course, free 
to choose a convenient set of p, u interpolating fields and work with these, provided one 
calculates contributions to S-matrix elements. Since, to Oivrid — m u ) Eqn. (2.2) remains 
valid when we replace p and to with and u/°\ the fields 



(0)c = 9p_y P 

rp " m 2V 



(0)c 9^yu ( 2A) 

satisfy(0|p^ c |p^ ^) = e M and (0\uj^ c \uj^) = e M , and hence serve as possible choices of in- 
terpolating fields for p^ and u/°). With this choice of interpolating fields (and not with 
others) one obtains 

IW) = ^A^(g 2 ) (2.5) 

where A^ pw (g 2 ) is the off-diagonal element of the vector meson propagator for the inter- 
polating field choice above. If one simultaneously evaluates, e.g. NNp, NNu vertex form 
factors using the same interpolating fields, one could, of course, piece the resulting vertices 
and propagators together to obtain contributions to NN scattering S-matrix elements which 
are independent of field choice. For a general choice of interpolating field, however, neither 



nor 11^ is proportional to A^. Given the existence of QCD sum rules and ChPT 
methods, which are rather efficient at handling current-current correlators and vector cur- 
rent vertex functions, the choice (2.4) for vector meson interpolating fields would appear 
to be a convenient and sensible one. With this choice, Eqn. (2.1) provides the basis of a 
spectral representation of A^ pul , but for other choices of the vector meson interpolating 
fields this is not the case. 

Note that the above discussion also clarifies one ongoing point of debate in the literature, 
namely that concerning the g 2 -dependence of the quantity 9 pul {q 2 ) appearing in Eqn. (1.3). 
Defining II(g 2 ) by 

iV(g 2 ) = (g»„ - Wg 2 )n(g 2 ) , (2.6) 

the absence of massless singularities implies that 11(0) = . This in turn implies, with 



A^(g 2 ) = -(g, u - Wg 2 )A^(g 2 ) , (2.7) 

A^(q 2 ) = 0, and hence 9 puJ (0) = 0. Since this is true for one choice of the vector meson 
interpolating fields, it is incumbent upon those advocating 

0p"(q 2 ) = 9 puJ (ml) (2.8) 

to explicitly demonstrate the existence of an interpolating field choice for the vector mesons 
for which Eqn. (2.8) is valid; the relation cannot be true in general. 

III. THE QCD SUM RULE ANALYSIS OF IV (Q 2 ) REVISITED 

With the above discussion in mind, let us turn to the sum rule analysis of the vector 
correlator, first briefly reviewing the treatment and results of Refs. [p~2|JT7|1 . The sum rule 



approach consists of writing an operator product expansion (OPE) representation for the 
correlator, valid in the region of validity of perturbative QCD, and a second, phenomeno- 
logical, representation in terms of hadronic parameters, and then Borel transforming both. 
The Borel transform serves to extend the ranges of validity of both representations and, in 
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addition, to (1) emphasize the operators of lowest dimension in the OPE representation and 
(2) give higher weight to the parameters of the lowest lying resonances in the phenomeno- 
logical representation. One then matches the transformed representations in order to make 
predictions for the relevant hadronic parameters. 

The OPE for the correlator of interest was performed long ago [jTjJ . Truncating the 
expansion at operators of dimension six, one finds that, defining IT(g 2 ) by 



n 



(3.1) 



one has 



n 

where Q 2 = —q 2 and 
c 

Cl 

c 2 
c 3 



(Q ) = ^ 



1 

12 



Cl 



c 2 , 2c 3 



OCEM 



16tt 3 



: (m 2 d - m 2 u ) 



2tt 2 
m d - m u 
m d + m u 
224 



7T [a,(qq)o] 



1 + 



OlEM 



7 



7 



m d + m u 



2 + 7/ Vm d - m u 



(3.2) 



(3.3) 



?a s (/i 2 ) 

with 7 = (dd) o / (uu) o — 1. Taking for the phenomenological representation (in the narrow 
resonance approximation) 



Imn phen (g 2 



7T 

12 



/ p % 2 - m 2 ) - /^(g 2 - m^) + f p ,5(q 2 - m 2 ,) - f u ,5(q 2 - mi,] 
a em 

1927T 2 ' 



+ 

(3.4) 



(where f p , f u , f p >, f^i may be thought of as the parameters to be determined from the sum 
rule analysis) one finds, upon Borel transformation and matching, 

1 



M 2 



/ p exp(-m 2 /M 2 ) - ^exp(-m 2 /M 2 ) + /„, exp(-m 2 /M 2 ) - ^,exp(-m 2 ,/M 2 



i67r3 exp(-WM) = c + — + — + M6 , 



(3.5) 



where M is the Borel mass. As pointed out in Ref. |12j , to 0(5m 2 , 5m' 2 ), where 5m 2 = 
m^j — m 2 , 5m' 2 = m 2 ^, — m 2 ,, Eqn. (3.5) can be rewritten in terms of the parameters £, /3, £' 
and /?', where 

5m 2 ( f p + f u \ 



m 4 \ 2 J 

(/» - fp) 

m 2 £ 
Sm' 2 ( fp' + U 



m' 4 V 2 



ft = l2 f ;' ] (3.6) 



with m 2 = (m^ + m 2 J )/2 and m' 2 = (m 2 ,/ + m 2 J ,)/2, as 

^ - ex P(--VM 2 ) + K - /?') exp(-m' 2 /M 2 ) 

+ 1^3 ex P (-,o/M )=co + W2 + W4 + j^. (3.7) 

If c _3 were precisely known, Eqn. (3.5) or Eqn. (3.7) could, in principle, be used to deter- 
mine the parameters £, (3, (3' . There are, however, some uncertainties in the values of 
the q, associated with the imprecision in our knowledge of the values of the four-quark con- 
densates and of the isospin-breaking ratio of the (uu)o and (dd)o condensates. The authors 
of Ref. [|l2j (which updates Ref. [0 ) consider a range of possibilities for these quantities, 
and also take for r = (m<2 — m u )/(md + m u ) the value r = 0.28, obtained from an analysis 
of pseudoscalar isomultiplet splittings |24j] employing Dashen's theorem |25] for the electro- 



magnetic contributions to these splittings. The last ingredient of the analysis of Ref. (12 
is the imposition of an external constraint on the hadronic parameter £, based on the ob- 
served interference in the p-u interference region in e + e~ — > 7r + 7r~. This constrained value, 
£ = 1.13 x 10 -3 , is based on (1) the assumed connection between the correlator and the 
propagator (presumably valid for the essentially on-shell value of the mixing, though not 
elsewhere) and (2) the assumption that direct -> 7T7T contributions to e + e — > ty + tt can 
be neglected (see Ref. [^BJ for a discussion of these issues). There appears to be no particu- 



larly good reason for the latter assumption, and, indeed, it would seem appropriate to allow 
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£ to be fit by the sum rule analysis as a test of this assumption (as will be done below), but 
let us follow the analysis of Ref. [I2[ for the moment. Using the sum rule, Eqn. (3.7), and 
imposing the constraint £ = 1.13 x 1CT 3 , as discussed above, the authors of Ref. |12| solve 
for /3, £' and f3' for four different input sets {q}. Using the expression (3.4) for ImIP' ieri (5 2 ) 
and the fact that Il(g 2 ) satisfies an unsubtracted dispersion relation, one may show that, to 
first order in 5m 2 and 5m' 2 , 

Ren(0) = -^[£(l-/?)+£'(l-/3')] . (3.8) 

Using the values of the parameters obtained in Ref. , one finds that the ratios of the 
contributions to Re 11(0) from the p'-u' region to those from the p-u region are 1.8, 0.8, 0.3 
and 0.8 for input sets I, II, III, IV, respectively. The failure of the results to be dominated 
by the nearby (p, u) singularities suggests that the phenomenological form employed for 
the spectral function may well be incomplete, either in missing low-lying contributions or in 
failing to include the effect of even more distant singularities. If we consider Eqns. (3.4) and 
(3.8) for a moment an interesting possibility becomes evident. If one had all isospin-breaking 
effects generated solely by p^-tu^ mixing, and if the physical vector mesons were a simple 
rotation of the isospin-pure basis (not in general true when the wavefunction renormalization 
matrix of the system is non- diagonal), we would have f p = f u for f p , f u as written in 
Eqn. (3.4). While the assumptions required to arrive at this conclusion are certainly not 
satisfied in general, this nonetheless indicates that there should be significant cancellation 
between the p and uj contributions to the correlator. Thus, a single isolated resonance, even 
with a coupling much smaller than that of the p or u, could in fact contribute significantly 
to U^,,. This suggests that the contribution to ImlLj,,, neglected in Ref. [0, may well be 



non- negligible. In fact we can make a rough estimate of the expected size of (where is 
defined by adding a contribution j^f^^q 2 — mi) to ImIT p?teri (g 2 ) in Eqn. (3.4)) as follows. 
<f) is known to be not quite pure ss. If, e.g., we take the Particle Data Group (PDG) P7| 
value for the octet-singlet mixing angle, 9 = 39° (quadratic fit), ~ — 5tu^°\ where 
is the pure ss state and 5 = .065 rad is the deviation of 9 from ideal mixing. The 
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contribution of the (f) pole term to IT^ due to mixing in the propagator should then be of 
order —5 times that associated with the u pole, i.e. ~ 0.065 / w ~ 0.065 f p . There will, 
of course, also, in general, be isospin-breaking contributions from direct couplings to the 
current vertices, not just from mixing in the propagator, but the above discussion shows 
that ~ (0.05 — 0.10) f PtU} should be a reasonable expectation. As we will see below, this 
(rather crude) estimate is indeed borne out by the sum rule analysis. 

Let us, therefore, add a term ^f<j,S(q 2 - m|) to lmll phen (q 2 ) on the RHS of Eqn. (3.4), 
and perform a reanalysis of that equation. We will follow Ref. Jl2[ in choosing the range of 
input values for the {q}, with, however, the following modifications. First, the small c\ term 
dropped in Ref. |12j will be retained, though, as pointed out there, it in fact has little effect on 



the final results. The numerical value is obtained by using (md+m u ) (lGeV) = 12.5±2.5MeV 
from Ref. and the updated value of r discussed below. The main modification to the 



input is in the parameter r. There is now considerable evidence that Dashen's theorem is 



significantly violated p9|-|31"|| , Refs. pC||3~Tf in particular suggesting that 



(m K+ - m K0 ) E M - 1-9 (m^+ - m^) expt (3.9) 

(where the factor 1.9 on the RHS of Eqn. (3.9) is absent in Dashen's theorem). Using 
Eqn. (3.9) in place of Dashen's theorem for the electromagnetic contribution to the kaon 
mass splitting produces a rescaling of r by 1.22. The resulting change in the Cj is essentially 
to rescale the values of Ci in Ref. |12| by this same factor. In assessing the effect of the 
uncertainties in the values of the {cj} for a given input set, the input errors on C2 have also 
been rescaled by this factor of 1.22. Finally, since the masses of all the resonances appearing 
above, including the p' and u/, are known, we may take these as input and use the sum rule 
to extract the isospin-breaking parameters, {fk}, where i = 1 • ■ - 5 correspond to p, uj, p', uj' 
and (J), respectively. Note that, in taking this approach, we are abandoning the constraint 
on £ employed in Ref. ||12|| . If the direct u>^ — > tt + it~ coupling is, indeed, negligible in 



e + e — > 7T + 7r , this will manifest itself by the value of £ resulting from the sum rule analysis 
being near 1.13 x 10~ 3 . 
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The analysis of the modified version of the sum rule, (3.5), proceeds as follows. First, 
from the terms of O(M ), cq = OEjvf/167T 3 . One may check that, as in Ref. |12j , the analysis 
is very insensitive to the value of the EM threshold parameter, sq. We will, therefore, quote 
all results below for the value, so = 1.8 GeV, employed in a number of the results quoted in 
Ref. |T2J . Second, again as in Ref. JT2J , we impose the local duality relation 



o 



oo 

ds ]mU phen (s) = 0{a EM ,m 2 g ) (3.10) 



(which is equivalent to matching the coefficients of the 0(1/M 2 ) terms in Eqn. (3.5)). With 
the index k — 1, • • • , 5 labelling p, lu, p', to' and 0, respectively, as above, this relation is 

£(-l) fc+ 7fc = Co*o + ci . (3.11) 

k 



(Note that the q tabulated in Ref. |T2 have had the appropriate factors of m 2 required 



to leave the remaining coefficient dimensionless factored out of them. Thus, e.g., C\ in 
Eqn. (3.11) is m 2 times that tabulated in Ref. [12| .) The remaining four relations required 
to obtain a solution for the five unknowns, {fk}, are obtained by acting on Eqn. (3.5) with 



v — ±j gnjjtfzyi for n — 1, • • • , 4. One may check that the results are not sensitive to using 
precisely the PDG values for the p' and u' masses. Indeed, shifting either mass by 50 MeV 
induces changes of < 4% in £, < 2.5% in /3, < 5% in (3' and < 20% in The resulting 
changes in the correlator itself are even smaller: e.g. 11(0) and Jpr(0) are changed by < 2% 
by the above mass shifts. 

In Table 1, the results of the modified sum rule analysis are displayed for the input sets I, 



III, IV of Ref. |T2[ , modified as described above. The errors shown in the table correspond 
to the uncertainties in the input parameters, c 2 and c 3 , (those quoted in Ref. [0 in the 
case of C3 and the rescaled version thereof in the case of c 2 ). The stability of the analysis 
is illustrated, for input set IV, in Figs. 1-5, which display the parameters £, (3, j3', 
as a function of the Borel mass, M, in the range 1-10 GeV (the choice of the first four 
parameters, rather than corresponding fk values, is made in order to facilitate comparison 
with Ref. Hm ). Set I generates results of comparable stability, while the results of set III 
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are even more stable than those of set IV. In all three cases a wide stability window exists 
in the Borel mass for all five output parameters. This stability window, moreover, occurs 
without the necessity of using unphysical values for the the average of the p' and u' masses. 
As noted previously in Ref. [12| , results for input set II are considerably less stable than 



for the other sets: in fact, no stability window exists anywhere in the range M — 1 and 
M = 10 GeV, apart from for the very lower edge of the error band for the magnitude of C3, 
for which values input set II is very close to the upper end of the corresponding error band 
for input set I. The instability of the analysis for input set II is illustrated (for the central 
values of c 2 and C3) in Fig. 6, where the parameter, fy, is plotted as a function of the Borel 
mass, M. As a result of this instability, results corresponding to input set II are not quoted 
in the table; for most of the input range (i.e. for larger values of the magnitude of C3) the 
input set appears, from the sum rule analysis, to be unphysical. 

A number of features are evident from the results of the above analysis. First, from Ta- 
ble 1, we see that the magnitude of £ differs significantly from that which would be expected 
from the analysis of e + e~ — > neglecting cu^ — > contributions, suggesting that 

the latter are, indeed, not negligible. It should be stressed that the errors quoted in the 
table correspond to varying c 2 and c 3 separately within the range of quoted errors, and 
taking the maximum variation of the resulting output. One can obtain even lower values 
of £, i.e. closer to that expected if one can indeed neglect cu^ — > ix + ix~ contributions to 
e + e~ — > 7r + 7r~, by letting c 2 lie at the bottom of its error band and, simultaneously, the 
magnitude of C3 lie at the top of its error band in set I. However, such a combination (which 
produces £ = 1.43 x 10~ 3 ) is quite unstable, the values of e.g., varying by more than 
20% between M = 3 and 5 GeV. A similar result, £ = 1.48 x 10~ 3 , can be obtained from 
set II for the central value of c 2 and the lower edge of the error band for the magnitude 
of c 3 , with comparable (~ 20% over the range M = 3 to 5 GeV) instability. All other 
portions of the set II error band are even more unstable. Thus it appears very clear that the 
value £ = 1.13 x 10~ 3 is excluded by the sum rule analysis. The second observation is that 
the inclusion of the pole term in the phenomenological representation of the correlator 
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rectifies the problem of the strength of the distant singularities. This can be seen from 
the relative size of £ and £' in Table 1, but is more evident in Table 2, where the output 
values for the parameters {fk} are tabulated, for the central values of the input parameters 
{cj}, for input sets I, III, IV. The ratios of to f u are 0.062, 0.068 and 0.066 for sets I, 
III and IV, respectively. This is in (better than should be expected) agreement with the 
rough estimate given above, confirming the physical plausibility of the solutions obtained. 
Moreover, f p i and f w > are now a factor of 40-60 smaller than f p and f w . The structure of 
the resulting contributions to the correlator near q 2 = is shown in Table 3, where the p, 
u), and also the p', u' contributions have been combined. Note that the individual p and 
u) contributions are a factor of ~ 13 larger than the contribution, but the cancellation 
between them is such that the <fi contribution is approximately twice as large as their sum. 
The p'-u' region contribution is then less than 10% of the <fi contribution. The more distant 
singularities, thus, have only a small effect, justifying, a posteriori, the neglect of yet more 
distant singularities in the phenomenological side of the sum rule analysis. Given that the 
results satisfy all the above tests for being physically sensible and stable, it appears that the 
resulting values for the correlator and its slope with respect to q 2 at q 2 = should be taken 
as good estimates, within the uncertainties resulting from those in the input parameters. 
The fact that, due to cancellation between the otherwise dominant p and uo contributions, 
the contribution is actually dominant, no doubt accounts for the unphysical behavior of 
the spectral distribution of the correlator obtained in the absence of the term. Note that, 



despite the significant changes in the fit, as compared to Ref. |T2[ , the slope of the correlator 
remains large in the present results. 

IV. THE CORRELATOR TO ONE-LOOP ORDER IN CHPT 

The starting point for the computation of the mixed-isospin correlator, ILj„(g 2 ), is the 



effective chiral Lagrangian of Ref. |2T 

C eff = ^fTr(D^D^) + i/ 2 Tr[ J B M(S + + ^[Tr^E^E*)] 2 
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+ L 2 Tr(D l j:D^)Tr(D^D v ^) + L 3 Tr(D^ D^ZD U ^ D U Z) 

+ L 4 Tr(^SD' i E t )Tr[2 J BoM(S + £ f )] + L 5 Tr[2 J B (MS + Y)M)D^D»Y] 

+ L 6 [Tr[25 M(S + £ f )]] 2 + L 7 [Tr[2fi Af (E - £ f )]] 2 

+ L 8 Tr[45 2 (MEMS + MS t MS t )] - iL 9 Tx[F R D^D v ^ + F^D^Y) ' D v Y\ 
+ LioTrp^^SF^] + H{Tx[F R F L ^ + F^F L ^} + # 2 Tr[45 2 M 2 ] . (4.1) 

In Eqn. (4.1), B is a mass scale related to the value of the quark condensate in the chiral 
limit, S = exp(zA ■ 7?//) (with A the usual SU(3) Gell-Mann matrices and 7? the octet of 
pseudoscalar (pseudo-) Goldstone boson fields), / is a dimensionful constant, equal to f n in 
leading order, M is the current quark mass matrix, and is the covariant derivative 

= 9 M £ - i(v„ + a M )S + zE(^ - a M ). (4.2) 

In Eqn. (4.1) the external pseudoscalar sources (which occur in the most general form of 
J~- e ff) have already been set to zero, and the external scalar source to 2B times the current 
quark mass matrix, since we are interested here only in the vector current correlator and, 
therefore, require only the external vector sources. For this same reason we may drop the 
external axial vector sources from the expression for the covariant derivative in Eqn. (4.2). 
The left and right external source field strength tensors, F^ R , then both reduce to Fz) R = 
dfj,v u — dyV^ — i[i> M , v u ], where v M = ^-v*, with the octet of external SU (3) vector fields. In 
principle, Eqn. (4.1) should be supplemented by terms involving Tr{F llu ) in order to treat 
the case at hand, since the current VJf contains both octet and singlet pieces. However, 
to one-loop order, the additional terms do not contribute to the isospin-mixed correlator 
(the correlator is identical to Uf u {q 2 )/3V3, with Uf u (q 2 ) = i J d*x(0\T(V*(x)V*(0))\0), to 
this order), so we will not explicitly display these terms. The unrenormalized higher order 
coefficients L 1; ■ • • , L w and Hi, H 2 appearing in Eqn. (4.1) contain divergent pieces which 
cancel those of the one-loop graphs involving vertices from the first two terms in £ e ff above, 
and also finite, renormalization-scale-dependent pieces, L\. Expressions for the divergent 



pieces of the Lj, Hi, relevant to one-loop calculations, may be found in Ref. |21 
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Contributions to the correlator resulting from Eqn. (4.1) are of two types, corresponding 
to the two types of contribution to the low-energy representation of the product of currents, 
V^Vy'- (1) those terms arising from the product of the low-energy representations of the 
individual currents, V£ and V" (obtained from the terms in C e ff linear in the v®, a = 
0, • • • , 8), and (2) contact terms (generated by the terms in C e ff quadratic in the v®). To 
leading order (i.e. keeping only the first two terms in £ e ff) the correlator vanishes. This is 
because it is isospin-breaking and the only isospin-breaking at leading order lies in the term 
involving the quark mass matrix, which does not contain the external vector sources, and 
hence does not contribute to the correlator in zero-loop graphs. The leading contributions 
to the correlator are, therefore, next-to-leading order in the usual chiral counting. As such, 
the contributions consist of one-loop contact and non-contact graphs (where the current 
vertices are obtained from the first term in the effective Lagrangian above) and meson-field- 
independent contact terms from the remainder of C e ff- These latter contributions, which 
would in general produce terms involving the Z£, may be easily shown to vanish for the case 
at hand. Thus only the contact and non-contact graphs mentioned above contribute. It is 
straightforward to demonstrate then that, to one-loop order, the Oim^ — m u ) expression for 
the correlator is 



log(m 2 K0 /m 2 K+ ) ( 4m| l\ - , / 4m|+ l\ - , 
48^ JK0{q) -[-^-3 JK+{q) 



(4.3) 



J P (q 2 ) = -/ dx log 1 - x(l - x)q 2 /m 2 P . (4.4) 



where 



I67H Jo 

For our purposes we will not need the general expression for J (which is quoted in Appendix 
A of Ref. [ETJ ), but only the behavior near q 2 = 0, which is given by 



1 q , 1 q_ 



Jp (q) = t^ti^y + t^i^t + ■■■ ■ ( 4 - 5 ) 

967r 2 m P 9607T 2 mj 



In Eqn. (4.3), Ttr^o K+ are the leading-order expressions for the kaon squared-masses, m 2 K o = 
Bq{itl s + rrid) and m 2 K+ = Bo(m s + m u ) and terms have been kept only to 0(irid — m u ). As 
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such, we must also expand all terms occuring there to the same order. Doing so, and making 



where fh 2 K is the average of the K + and K° squared masses. Thus, 1211(0) = {m 2 K0 — 
m 2 K+ ) l '487r 2 m 2 K , where the kaon mass difference is that due to the strong isospin-breaking, 
i.e., with the electromagnetic contribution removed. Using Eqn. (3.9) for the electromagnetic 
contribution, we find that the RHS of this expression is 5.5 x 1CT 5 , to be compared with the 
results of the sum rule analysis, ~ 1 x 1CT 3 . The one-loop ChPT result is a factor of ~ 20 
smaller than the sum rule result. 

The discrepancy between the one-loop ChPT and sum rule analyses for the correlator 
near q 2 = should actually not come as a complete surprise. Indeed, when the leading- 
order contribution to a physical quantity (order 2 in the chiral expansion) vanishes, as it 
does here, one has no obvious scale to use in judging whether or not the next-to-leading- 
order contribution obtained is abnormally small, i.e., whether or not the resulting one- 
loop expression is likely to represent a well-converged approximation to the whole chiral 
expansion. In fact, the structure of the expression, (4.3), above for the correlator, II(g 2 ), is 
such as to suggest that it is unlikely to be well-converged. The reason for this statement 
is that Eqn. (4.3) is independendent of the low-energy constants (LEC's), L£, and results 
purely from one-loop graphs involving internal kaon loops. Such loops, for the non-contact 
graphs, are well-known to be suppressed in size (the coefficient of q 2 in the leading term 
of Jk in Eqn. (4.5), e.g., is a factor of m 2 /m|- smaller than for the corresponding 7r loop 
integral, J^) and, moreover, in the case at hand, i.e. the correlator IT(g 2 ), those terms in 
which this suppression would be lifted by the presence of the m 2 K /q 2 factor in the coefficient 
multiplying Jft-(g 2 ) cancel, since the expression for the correlator involves the difference 
of the K + and K° loop contributions. The correlator, of course, has a cut beginning at 
q 2 = 4m 2 , associated with nn intermediate states, but such intermediate states do not enter 



the expansion of Eqn. (4.5) for the loop integrals Jr°,k+, we obtain, for the behavior of the 
correlator in the vicinity of q 2 = 0, 




(4.6) 
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until two-loop order in the chiral expansion. Since the relevant 7r loop integral is intrinsically 
much larger than its kaonic counterpart, it is likely that the two- loop contributions will not 
be negligible, despite being higher order in the chiral expansion. 

Other examples of slow convergence of the chiral series when the leading contribution 
vanishes and the next-to-leading order contribution results purely from loop graphs (i.e. is 
independent of the fourth-order LEC's, L£) are, in fact, already known. One is the process 
77 — > 7r°7r°, whose amplitude, to one-loop order, receives contributions only from loop 
graphs (though in this case, loop graphs with internal 7r lines). The one-loop expression 
|32| , |33H deviates from the experimental amplitude |34| even very close to threshold, and one 
finds that extending the calculation to two-loop order (sixth order in the chiral expansion) 
produces corrections to the one-loop result of order 30% |H| , which corrections bring 
the amplitude into agreement with experiment. Even more closely similar to the case at 
hand is the process r\ — > 7r°77. The one-loop amplitude again has no leading term and no 
contributions from the fourth order LEC's, but here, although there are 7r loop contributions, 
these contributions are suppressed by a factor (m^ — m u ). The K loop contributions are 
naturally small, as noted above. The result is that the one-loop prediction for the partial 



rate is a factor of ~ 170 smaller than observed experimentally 

It is worth considering the process rj — > 7r°77 in somewhat more detail since, not only 
does it closely parallel the case at hand, but the physical origin of the smallness of the 
one-loop result for this process is well-understood. The source of the problem lies in the 
fact that the dominant contribution to the amplitude is known to be due to vector meson 
exchange . As is well-known p],^,f|(| , it is possible to make standard field choice 



for the various meson resonances and write an effective chiral Lagrangian which includes 
both these resonances and the octet of pseudoscalar (pseudo-) Goldstone bosons. One may 
then integrate out the (heavy) resonance fields to obtain an effective Lagrangian of the 
form C e ff for the pseudoscalars alone. The resonance contributions to the LEC's are then 
determined by the coupling parameters of the original, extended Lagrangian (which are 
fixed by experimental data). The effect of the resonances (for the initial field choices used) 
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then lies solely in their contributions to the L£ fl3~9"|,[40"l| . Those L\ to which the vector 



and axial-vector mesons can contribute (Z£, % = 1,2,3,9,10) are known to be essentially 



saturated by these contributions |p9| , pJ|| . Thus, the absence of any Independence in the 
one-loop amplitude implies the absence of the effects of vector meson exchange (for the 
given initial choice of vector meson interpolating fields) and if, as seems to be the case for 
77 — > 7r°77, vector meson exchange is the dominant contribution, the one-loop amplitude 
can be expected to be a poor representation of the full chiral series. The vector meson 
contributions, in this case, first appear as tree-level contributions arising from the 0(p 6 ) 
part of the effective Lagrangian, not included in Eqn. (4.1) above (the general form of the 
0(p 6 ) part of the effective Lagrangian is given in Ref. |Q ). Thus, only by including these 
contributions (which requires, for consistency, a full two-loop-order calculation) can one 
hope to obtain a well-converged approximation to the full chiral series for the amplitude. 

The 7] — > 7r°77 discussion above can obviously be transferred directly to the case of the 
mixed-isospin vector current correlator under consideration here. We expect significant con- 
tributions to the correlator from the vector meson resonances and, for a particular choice of 
vector meson interpolating fields, these contributions are completely absent from the one- 
loop result. As a consequence, we can expect significant contributions from the tree-level 
0(p 6 ) terms in which such contributions reside. As already discussed, the 0(p 6 ) loop contri- 
butions (arising from two-loop graphs with lowest-order vertices and one-loop graphs with 
a single 0(p 4 ) vertex) may also be significant. A two-loop calculation is, therefore, almost 
certainly required in order to obtain convergence of the chiral series for the mixed-isospin 
correlator. Similar statements hold for the related correlator, U 3 ^(q 2 ), for which work on the 
two-loop calculation is in progress [[42] . Note that, in the latter case, only a single combina- 



tion of the 0(p ) LEC's enters the two-loop result. This combination, which in the notation 
of Ref. jKJ (where the analogous IT^ and II®® correlators are computed to two-loop order), 
is written Q°(fi) — 3Lg ^(/i) — 3Li ^(/i), with fi the renormalization scale, is in principle 
obtainable from experimental data using the chiral sum rules of Ref. [^3| (Eqns. (97) and 
(98) therein). In the case at hand, one further 0(p 4 ) LEC and one further 0(p 6 ) LEC will 
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be present at the two-loop level, but the sum rule analysis above, in combination with a 
full two-loop evaluation of U(q 2 ), would provide a useful constraint on these parameters, 
albeit it with the ~ 20 — 30% errors displayed in Table III and associated with the un- 
certainties in the values of the input parameters {q} which determine the correlator near 
q 2 = 0. A similar sum rule analysis of the correlator nj^ would constrain the combination 
Q°(fi) — 3^^(12) — 3Liq (A*) > mentioned above. 

It is interesting to note that the relation between the sum rule and ChPT results for the 
mixed-isospin vector correlator is effectively the reverse of what occurs in the mixed-isospin 



axial correlator case. In the latter case, the ChPT |16j and sum rule |L3| results for the 
value of that piece of the correlator proportional to q^q u at q 2 = are comparable, but 
the ChPT result for the slope of the correlator with q 2 is more than an order of magnitude 
larger than that obtained from a sum rule analysis analogous to that employed above for 



the vector correlator case [16[ . The source of the discrepancy, in the axial correlator case, is 
that the sum rule result for the slope has the incorrect chiral behavior, being in fact missing 
its leading contribution in the chiral expansion. This problem with the sum rule treatment 
is easily exposed using chiral methods, but is completely non-obvious without them. In the 
present case, since we do not know what portion of the vector meson masses survives in the 
chiral limit, we cannot make as precise statements about the required chiral behavior of the 
vector correlator. There is, however, no obvious problem with the form of the sum rule result 
above. The sum rule result, moreover, provides clear evidence to indicate that the chiral 
series for the vector correlator is indeed, as suggested by analogy to the known behavior 
of the i] — > 7r°77 process, slowly converging. The sum rule result, in this case, should also 
provide useful input for the two-loop analysis in ChPT. The two examples clearly indicate 
the advantages of applying both methods, within their common range of validity, in any 
given physical process. 
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V. SUMMARY OF RESULTS 



The basic results of the paper are as follows. We have demonstrated that (1) in making 
a sum rule analysis of the mixed-isospin vector current correlator, it is necessary to include 
the 4> pole term in the phenomenological form of the representation of the correlator, and 
that, when one does so, the spectral structure of the correlator becomes physically sensible; 
(2) the expression for the correlator away from q 2 = has no general interpretation 
as the off-diagonal element of a vector-meson propagator except for a particular vector 
meson interpolating field choice; (3) the freedom of field redefinition shows that the isospin- 
breaking factor 9 pu '(q 2 ), which occurs in the numerator of the expression, (1.3), for the 
off-diagonal element of the vector meson propagator, cannot, in general, be taken to be 
independent of q 2 ; (4) the behavior of the correlator near q 2 = suggests that the direct 
o/°) — > 7r + 7r~ contribution to e + e~ — > 7r + 7r~ is not negligible in the p-u interference region; 
(5) the discrepancy between the behavior of the correlator near q 2 = as obtained from 
the sum rule analysis and from ChPT to one-loop indicates a slow convergence of the chiral 
series for the correlator and, in consequence, the necessity of a two-loop calculation of this 
quantity in ChPT. The sum rule result for the correlator near q 2 = can then be used, in 
such a calculation, to constrain the 0(p 6 ) LEC's of ChPT. 
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FIGURES 
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FIG. 1. Dependence of £ on the Borel mass, M, for modified input set IV. 




FIG. 2. Dependence of (3 on the Borel mass, M, for modified input set IV. 
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FIG. 3. Dependence of £' on the Borel mass, M, for modified input set IV. 
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FIG. 4. Dependence of (3' on the Borel mass, M, for modified input set IV. 
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FIG. 5. Dependence of on the Borel mass, M, for modified input set IV. 




FIG. 6. Dependence of on the Borel mass, M, for modified input set II. 
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TABLES 





TABLE I. 


Sum rule fit for the parameters £, 


P, P' and / 




Input 


£ x 10 3 


P 


£' x 10 5 


P 


U x io 3 


Set I 
Set III 
Set IV 


2.18±0.39 
3.10±0.39 
2.59±0.39 


1.49L0.06 
1.62L0.02 
1.55L0.04 


-2.63±0.79 
-4.57±0.69 
-3.47L0.61 


-5.84±0.12 
-5.72L0.01 
-5.78±0.04 


2.30±0.52 
3.57±0.52 


TABLE II. Sum rule fit for the isospin-breaking parameters {/&}. Values are 
central values of the input parameters {cT. The units are GeV 2 . 


quoted for the 


Input 


f P x 10 2 


f u x 10 2 


U x 103 


fp> x 10 4 


fw> x 10 4 


Set I 
Set III 
Set IV 


3.53 
5.00 
4.18 


3.73 
5.30 
4.42 


2.30 
3.57 
2.86 


5.34 
9.32 
7.06 


8.45 
14.6 
11.1 


TABLE III. Behavior of the correlator near q 2 = 0. Contributions to 1211(0) from the p-u, 
4> and p'-u)' regions are quoted for central values of the input parameters {q} for each input set, 
while the effect of the uncertainties in these values is displayed explicitly for 1211(0) and 12^(0). 
All entries are in units of 10~ 3 , except for 12^(0), which is in units of 10~ 3 GeV -2 . 


Input 


p-LO 


<P 


p'-u 1 


1211(0) 




Set I 
Set III 
Set IV 


-1.06 
-1.92 
-1.43 


2.21 
3.44 
2.75 


-0.18 
-0.31 
-0.24 


0.96±0.14 
1.22L0.14 
1.08L0.14 


3.88±0.61 
5.10±0.60 
4.43±0.61 
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